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We study fidelity decay in classically chaotic microwave billiards for a local, pistonlike boundary 
perturbation. We experimentally verify a predicted nonmonotonic crossover from the Fermi golden 
rule to the escape-rate regime of the Loschmidt echo decay with increasing local boundary pertur- 
bation. In particular, we observe pronounced oscillations of the decay rate as a function of the 
piston position which quantitatively agree with corresponding theoretical results based on a refined 
semiclassical approach for local boundary perturbations. 

PACS numbers: 05.45.Mt, 03.65.Sq 



I. INTRODUCTION 

The stability of quantum time evolution measured by 
the overlap between time-evolved perturbed and unper- 
turbed states, as suggested by Peres [T| , has been studied 
from various viewpoints and under different names. In 
the field of quantum information this overlap is called 
"fidelity" Q and plays an important role for quantify- 
ing the susceptibility of quantum dynamics to environ- 
mental or other external perturbations. In semiclassical 
quantum and wave mechanics, alternatively, the overlap 
of an initial state with the state reached after succes- 
sive forward and backward time propagation, governed 
by the unperturbed and perturbed Hamiltonian, is often 
termed "Loschmidt echo" (LE) Q, especially for Hamil- 
tonians associated with complex, in particular, classically 
chaotic dynamics. This terminology refers to the notion 
of echoes from momenta reversal in a Hamiltonian system 
considered by Loschmidt [jj in the 19th century. 

For chaotic systems the LE has been predicted to ex- 
hibit different decay characteristics [5[ depending on the 
form and strength of the perturbation. One distinguishes 
roughly three prominent LE decay regimes: the per- 
turbative Gaussian [1,0, the Fermi-gplden-rule (FGR) 
JEHU, and the Lyapunov regimes d, (for reviews see 
Refs. [l3, EH)- The various perturbations considered 
have in common that they act "globally" on the system, 
i.e., already a moderate perturbation strength can cause 
a considerable rearrangement of the spectrum and eigen- 
functions. Correspondingly, in a semiclassical picture, 
a global perturbation affects all trajectories of the sys- 
tem, and hence all of them are responsible for the decay 
of the LE. The corresponding, original semiclassical ap- 
proach to the LE [j| 0] , which was recently generalized 
beyond the so-called diagonal approximation [l2j |. was 
extended in Ref. (l3j to strong local perturbations in co- 
ordinate space. This combined analytical and numerical 
study revealed for a billiard with a local boundary de- 



formation, much larger than the de Broglie wavelength, 
a novel LE decay law exp(— 2^t), where 7 is the classi- 
cal "escape rate" from the related open billiard. This 
approach was refined and generalized to weak perturba- 
tions in Ref. [l4| predicting a nonmonotonic crossover 
from the FGR to the escape-rate regime with increasing 
perturbation. For the case of a pistonlike boundary per- 
turbation the LE decay rate is expected to show distinct 
oscillations as a function of the perturbation strength, 
i.e., piston position. While this nonmonotonic crossover 
has been numerically confirmed for maps |l5j , quantum 
wave packet simulations for billiards requiring more ex- 
pensive numerics could only reveal precursors of this be- 
havior [l4| , calling for an experimental verification of the 
oscillations. 

For a global perturbation the fidelity decay was studied 
in a microwave billiard with classically chaotic dynamics 
by shifting a billiard wall [l6| . Using the concept of scat- 
tering fidelity [TtJ the predicted fidelity decay from the 
perturbative to the FGR regime was verified experimen- 
tally. 

A theoretical and experimental investigation of fidelity 
decay for another type of "local" perturbation in the 
perturbative regime, where the eigenstates are not sig- 
nificantly modified by the perturbation, has been done 
in [18[. On the experimental side a small scatterer was 
shifted inside the microwave billiard in a two-dimensional 
array of pointlike scatterers. Using the random plane 
wave conjecture, an algebraic decay 1/t was predicted 
theoretically and confirmed experimentally. Another 
type of fidelity decay caused by local perturbations has 
been studied in Ref. [lj|, where the coupling to an at- 
tached antenna was varied. In the present paper we use a 
microwave billiard with a piston attached to address the 
predicted nonmonotonic features in the fidelity decay. 

The paper is organized as follows. In Sec. [II] we present 
the experimental setup and introduce the scattering fi- 
delity. In Sec. Mil we briefly summarize the semiclassical 
results for the LE decay with local boundary perturba- 
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FIG. 1. (Color online) Geometry of the chaotic Sinai-shaped 
billiard (length of 472 mm, width of 200 mm, and a quarter- 
circle of radius 70 mm) with a variable pistonlike local bound- 
ary deformation. The piston position can be changed from a 
displacement h — 45 mm to h = mm for four different piston 
widths w — 20, 40, 70, and 98 mm. At position a the mea- 
suring antenna is introduced. The additional elements were 
inserted to perform ensemble averages (rotatable ellipse) and 
to reduce the influence of bouncing balls. 



tions [14| and derive (in the Appendix) an extension of 
the expression for the effective decay rate for the case of 
a pistonlike boundary perturbation as used in the exper- 
iment. We then present in Sec. IIVI our results for the 
experimentally determined scattering fidelity decay and 
compare them with the theoretical predictions for the 
corresponding Loschmidt echo decay. Our main findings 
are then summarized in Sec. [V] 



II. EXPERIMENT 

Microwave experiments with fiat cavities have become 
a well-known paradigm in the field of quantum chaos pol ] . 
In microwave billiards we can measure scattering matrix 
elements S a b(y) and S' ab (y) for unperturbed and per- 
turbed systems, independently, in frequency space. The 
scattering fidelity amplitude is defined in terms of their 
Fourier transforms, S a b an d S'„ h (upon choosing an ap- 
propriate frequency window) [17] : 



/«&(*) = 



(s ab (t)s': b (t)) 



(Sa b (t)s* ab (t)}(S' ab (t)S'; b (t)} 



The scattering fidelity itself is 

F(t) = \f ab (t)\ 2 . 



(1) 



(2) 



For chaotic systems and weak coupling of the measuring 
antenna the scattering fidelity approaches the ordinary 
fidelity 03 • 

In the present experiment we chose a resonator with 
a height of 8 mm which can be considered as two- 
dimensional for frequencies below 18 GHz. The setup, as 
illustrated in Fig. [TJ is based on a quarter Sinai-shaped 
billiard. Additional elements were inserted into the bil- 
liard to reduce the influence of bouncing-ball resonances. 
The classical dynamics for the chosen geometry of the 



billiard is chaotic. The straight left boundary of the un- 
perturbed billiard was deformed at a certain position by 
inserting pistons of four different widths w. The horizon- 
tal piston position can be changed in steps of 0.5 mm via a 
step motor from a displacement h = 45 mm to h — mm. 
At position a an antenna is fixed and connected to an Ag- 
ilent 8720ES vector network analyzer (VNA), which was 
used for measurements in a frequency range from 2 to 
18 GHz with a resolution of 0.1 MHz. We measured the 
reflection ^-matrix element S aa for four piston widths 
and all displacements h realizing 18 different positions of 
a rotating ellipse (see Fig. [1]) to perform ensemble aver- 
ages. The unperturbed system is defined as the one with 
the straight wall, corresponding to h = mm. 



III. THEORY 



The LE 



M(t) = m 



iH't/h -iHt/h\ 



(3) 



is defined as the overlap of an initial state \<p) evolved 
in time t under a Hamiltonian H with that evolved un- 
der a perturbed Hamiltonian H'. Within a semiclassical 
approach this quantity was studied in Refs. [HGJ] for lo- 
cal perturbations in chaotic systems. There it was shown 
that the LE is approximately [13] 



M(t) 



with the effective decay rate k given by 
k = 2 (l - Re(e 27r ™/ A )) 



(4) 



(5) 



Here u, called the deformation function, equals the length 
difference, induced by the local boundary perturbation, 
between the perturbed trajectory and the unperturbed 
one. A denotes the de Broglie wavelength. For the case 
of a pistonlike boundary deformation with piston width 
w and displacement h of the piston, as it is realized in 
our experiment, we find, in generalization of the results 
of Ref . [HI (see the Appendix) , 



w ^— ' 

k=0 



dx d6 cos 9 cos 



4-7T 

— (h cos 9 + kw sin 0) 
A 

(6) 

with the integration domains f^/c+i over incident posi- 
tions x and momentum directions 9 defined in Eq. (|12|l. 
In the limit h <C w, Eq. ^ reduces to [3] 



K = 7rHi(47i7i/A) 



(7) 



with Hi being the Struve ii-function of first order. 

Furthermore, M (t) in Eq. (j4|) depends on 7, which is 
the classical escape rate of the corresponding open cavity 
if the piston is removed. It is given by 



7 : 



Pa 
mid 



(8) 
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FIG. 2. (Color online) Measured scattering fidelity decay 
F(t), Eq. ([2]), (solid lines with symbols) for three different 
piston displacements hi = 1 mm (blue triangles) , hi = 5 mm 
(green circles), /13 = 10mm (red squares), for a frequency 
range 17 — 18 GHz corresponding to a mean de Broglie wave- 
length A w 17 mm. The dashed lines show the corresponding 
semiclassical prediction, Eq. Q, for the LE decay, with re cho- 
sen as free parameter: Ki — 0.26, K2 = 2.78, and K3 = 1.09, 
respectively. The time is given in units of the dwell time I/7, 
with 7 determined from experimental parameters via Eq. © 
with w — 40 mm. 



FIG. 3. (Color online) Decay rate function of piston 

displacement h for a piston of width w — 40 mm in a frequency 
range 17 — 18 GHz corresponding to a mean de Broglie wave- 
length A ~ 17mm. The asterisks represent the data points ob- 
tained from fitting the decay exponent of the measured scat- 
tering fidelity. The three cases discussed in Fig. [5] are marked 
by correspondingly colored symbols. The dashed curve shows 
the theoretical approximation (valid for h -C w), and the 
solid curve is a result of the numerical evaluation of the full 
semiclassical expression @. 



IV. RESULTS AND DISCUSSION 



for particles with momentum po and mass m, and for 
the average dwell length Id of paths in the related open 
chaotic billiard. In billiards with openings (deformation 
widths) w much smaller than the perimeter one can ap- 
proximate Id ~ ttA/w with A the area of the correspond- 
ing closed billiard. Further we will set po/m = c, where 
c is, in the case of the microwave billiard the speed of 
light. 

In semiclassical theory of the Loschmidt echo the per- 
turbation strength is a measure of the action change in- 
troduced by the perturbation. Thus, for perturbations 
caused by pistonlike boundary deformations the piston 
displacement h serves as the measure of the perturba- 
tion strength. As shown in Ref. [l4| it is convenient to 
define a dimensionless quantity \ = 27r-y/ (u 2 ) / A, where 
(u 2 ) = 8h 2 /3 for a pistontype deformation, as the pertur- 
bation strength. Then, based on Eqs. (Q]and[5]), different 
decay regimes of the LE can be identified as follows. For 
weak local perturbations, x < 1, one has M(t) » e~ x 7 * 
characterizing the FGR regime. Strong local perturba- 
tions, x ^* 1) l ea d to M{t) ~ e~ 2lt corresponding to the 
escape rate regime. In the following section we will use 
Eq. together with the refined expression ^ for the 
decay rate re, for a comparison with the experimentally 
determined scattering fidelity ©. 



In this section we present our measurements of the 
scattering fidelity decay for the pistonlike boundary per- 
turbation and compare them with the theoretical predic- 
tions (01)- © for LE decay for this specific type of per- 
turbation. We start with a piston of width w = 40 mm. 
In Fig. M the scattering fidelity F(t), Eq. ©, is plotted 
for three different piston displacements h acting as per- 
turbation to the system (symbols and solid lines). Ad- 
ditionally the corresponding semiclassical predictions for 
the Loschmidt decay according to Eq. (Q} are plotted 
(dashed lines) with re used as a fitting parameter, while 
7 was obtained from the geometry. The experimental 
fidelity decay shows good agreement with the expected 
exponential law beyond a certain time, which passes un- 
til the perturbation is "seen" during the measuring pro- 
cess. Upon increasing the displacements h, illustrated 
in Fig. [5] by the successive triangle (blue), circle (green), 
and squares (red) traces, the corresponding LE decay ex- 
ponent re exhibits a nonmonotonic behavior. 

For a more detailed investigation of this dependence of 
the Loschmidt decay exponent re on the displacements h 
of the piston, re is compared to the corresponding theo- 
retical predictions in Fig. [3] The data points re(ft), ob- 
tained from fitting to the experimental fidelity decay re- 
sults as in Fig. (5J are shown by asterisks. The three 
cases discussed in Fig. [2] are marked by correspondingly 
colored symbols at h\ =1, 5, and 10 mm. The figure 
shows an oscillating behavior of the fidelity exponent, 
and thereby of the fidelity decay at fixed time, as al- 
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FIG. 4. Decay rate K, ELS db function of x = \/8/3h 2ty/X. The 
fuzzy trace depicts the overlayed experimental data and the 
dashed curve the theoretical prediction (0). 



ready predicted and referred to as Fabry-Perot-type in- 
terferences between perturbed and unperturbed paths in 
14j . In Fig. |3] the semiclassical results for n(h) 



Ref. 

resulting from the numerical evaluation of the expres- 
sion © and the approximation ([7]) (for h <C w) are de- 
picted as solid and dashed curves, respectively. Already 
the dashed curve exhibits qualitatively good agreement 
with the measurement, though there is a mismatch in 
the amplitudes: While this approximative theoretical re- 
sult shows a monotonic decay of the maximum ampli- 
tude with h, the experimentally observed peaks of n(h) 
do not show this simple structure. However, the refined 
semiclassical prediction © (solid line) reflects the exper- 
imentally found irregular oscillation of amplitudes much 
more convincingly, showing reasonable agreement. In 
particular, for h « w = 40 mm (square shape of the 
pistonlike deformation) the experimental results show a 
particularly pronounced amplitude which is met by the 
solid line. In this range, which is beyond the range of 
validity of Eq. ([7]), the expression ^ constitutes a clear 
improvement. 

Furthermore, in Fig. @] we present on the same plot 
experimental data for k versus \ curves with the mean 
de Broglie wavelength A in the frequency range A < 2w, 
while w > 40 mm and h < w. We find very convincing 
agreement with the theoretical prediction (UJ). 

Finally, we demonstrate that the agreement between 
the experimental and theoretical curves can be shaken by 
pushing the experimental conditions too far beyond the 
main limit of the semiclassical theory, A <C w. Figure [5] 
shows the decay rate k for a piston width w = 20 mm 
which is of the order of A. As expected, the agreement 
between theory and experiment is not as good as that for 
the w = 40 mm case; in particular the experimental data 
points (dashed line) oscillate with a period that differs 
from the theoretical one. Experimentally, we again find 
a more pronounced amplitude around h ~ w, which is 



FIG. 5. Decay rate k as a function of the displacement h 
for a thin piston of width w — 20 mm in a frequency range 
17 — 18 GHz corresponding to a mean de Broglie wavelength 
A ~ 17 mm. The asterisks show the data points extracted 
again from the fit exponent of the exponential decay of the 
observed scattering fidelity. The dashed and solid curves show 
the theoretical predictions based on Eq. ([7]) and the numerical 
evaluation of Eq. ©. 



again described more convincingly by the full (solid line) 
than the approximative theoretical expression. However, 
the fact that the experimental parameters are beyond 
the regime of validity of the semiclassical theory does 
not allow for a further reasonable comparison between 
experiment and theory. 



V. CONCLUSIONS 

In this work we presented an experimental verification 
of the recent semiclassical predictions for fidelity decay 
arising from a local perturbation of a chaotic quantum 
system. In particular, we could confirm that the rate 
governing exponential fidelity decay exhibits oscillations 
as a function of the perturbation strength. The ob- 
served nonmonotonic behavior implies that for certain 
ranges of the perturbation strength the fidelity decay 
becomes weaker (for fixed time) with increasing pertur- 
bation strength. While the original semiclassical treat- 
ment [blj for a pistontype local boundary deformation 
was based on the assumption of a small piston displace- 
ment, the present microwave setting required a general- 
ization of the semiclassical approach to arbitrary ratios 
between piston displacement and width, which we per- 
formed by deriving an expression for the decay exponent 
in terms of a quadrature. We find quantitative agreement 
between the measurements and this refined semiclassical 
theory despite the fact that the microwave billiard does 
not really satisfy the underlying semiclassical assump- 
tion, namely, that the extent of the local perturbation, 
here the piston width w 7 should be much larger than the 
de Broglie wave length A. An improved semiclassical ap- 
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proach for local perturbations of size w < A would require 
one to treat semiclassical contributions due to diffractive 
trajectories properly, which is left for future research. 

On the experimental side, there remains the challenge 
to observe fidelity decay in the escape rate regime (for 
strong perturbations) characterized by a perturbation- 
independent fidelity decay rate. Naturally, this regime is 
difficult to access since the expected signals are tiny. 
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APPENDIX: SEMICLASSICAL THEORY OF 
EXPONENTIAL FIDELITY DECAY FOR 
ARBITRARY w AND h 

In view of the above experiments we extend the semi- 
classical theory for the decay of the fidelity due to pis- 
tonlike boundary deformations, presented in Ref. [lij for 
the limit h w. There it was shown that the rate of 
exponential fidelity decay is given by Eq. ([5]) with 

(e 2 ™/ A )= — / d9cos9e 27rm( - x ^/ x . (9) 
Jo w J 

Here x and 9 denote the incident position and angle, 
respectively [see Fig. [^a)]. The deformation function 
u(x,9) equals the length difference between the per- 
turbed and unperturbed trajectory of a correlated tra- 
jectory pair. A pair made up of an unperturbed and 
perturbed trajectory is considered correlated if the two 



FIG. 7. (Color online) "Unfolded" representation of corre- 
lated trajectory pairs belonging to sets fii (a), O3 (b), and 
ft n (c) with an odd integer n. 

trajectories exit the perturbation region with the same 
momentum direction 

QJ] (see Fig. El). 
We first note that only perturbed trajectories with an 
odd number of reflections may exit the perturbation re- 
gion with the same momentum direction as the unper- 
turbed trajectory and, therefore, contribute to the fi- 
delity. We denote by £l n a set of all correlated trajectory 
pairs where the perturbed trajectory exhibits n reflec- 
tions. The panels (a)-(c) in Fig. [5] show representative 
trajectory pairs belonging to the sets fii, and Q5, 
respectively. Equation (j9]) can then be written as 

00 „ 

(e 27rm/A ) = -V / dxd9cos9e iu/x , (10) 

where the double integral in the fcth summand runs over 
a region in the (x,9) plane that defines the set f^fc+i- 

In order to calculate the deformation function u(x,9) 
for a trajectory pair from the set fi ra we "unfold" the per- 
turbation rectangle (by "gluing" mirror copies of the rect- 
angle along the reflection sides) making the perturbed 
trajectory become a straight line (see Fig. [7]). Thereby, 
Fig. [7{a) is the "unfolded" version of Fig. [^a), Fig.^b) 
corresponds to Fig. EJb), and Fig. [7jc) represents a tra- 
jectory pair belonging to Q n . In this representation, the 
deformation function u equals the length of the interval 
AB, u = AB. Here A is a point of incidence, while the 
point B belongs to the perturbed trajectory and is spec- 
ified by requiring the angle ABA' to be tt/2, where A' 
represents the incident point A in the "exit" copy of the 
perturbation rectangle (see Fig. [?])• Then, a geometrical 
calculation yields 

u(x,6) =2/icos0+(n-l)wsin0, (II) 

where (x, 9) £ fl n . 

We now give a precise definition of the region fi n in 
the (x, 9) plane. As evident from Fig. Efc), a trajectory 
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FIG. 8. Schematic representation of the regions fii, Q3, and 
Qs,, see Eq. (|12[1 . Further regions, Q,2k+i with k > 3, con- 
tributing to the sum on the right hand side of Eq. (|6} are not 
shown in the figure; they cluster as narrow stripes "to the 
right" of SI5 and approach 9 = n/2 in the limit k — > 00. 



pair belongs to fl n if and only if the x coordinate of the 



exit point of the perturbed trajectory (in the unfolded 
picture) lies between (n — l)w and nw. This yields 

fl n = {(x,6) : xe(0,w), 8e (0,vr/2), 

x + 2htaxi9 e (ra- l,n)wj . (12) 

Figure [5] schematically shows the first three sets, fii, Q3, 
and O5, contributing to the sum in Eq. (p~0|) . 

Combining Eqs. ([5]), ([10]), (JTT]), and ([12]) we arrive at 
Eq. ([5]), which completes the derivation. The sum and 
integrals in Eq. are then computed numerically by 
means of Monte Carlo sampling. 

As a final remark, we note that in the limit of h <C w 
the sum on the right hand side of Eq. © is dominated by 
the k = term. Then the integration region tti can be 
approximately extended to the rectangle i£ (0, w), 9 £ 
(0, 7r/2), since the contribution for large angles, 9 close 
to 7r/2, is suppressed by the cos£? term in the integrand. 
This approximation leads to Eq. ([7]) (see also Ref. (3)- 
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